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Abstract 



We report on a new class of critical states in the energy spectrum of 
general Fibonacci systems. By introducing a transfer matrix renormaliza- 
tion technique, we prove that the charge distribution of these states spreads 
over the whole system, showing transport properties characteristic of elec- 
tronic extended states. Our analytical method is a first step to find out the 
link between the spatial structure of these critical wave functions and the 
quasiperiodic order of the underlying lattice. 
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1 



Typeset using REVI^X 



The notion of critical wave function (CWF) has evolved continuously since its intro- 
duction in the study of aperiodic systems^, leading to a somewhat confusing situation. For 
instance, references to self-similar, chaotic, quasiperiodic, lattice-like or quasilocalized CWFs 
can be found in the literature depending on the different criteria adopted to characterize 
themfM. Generally speaking, CWFs exhibit a rather involved oscillatory behavior, display- 
ing strong spatial fluctuations which show distinctive self-similar features in some instances. 
As a consequence, the notion of an envelope function, which has been most fruitful in the 
study of both extended and localized states, is mathematically ill-defined in the case of 
CWFs, and other approaches are required to properly describe them and to understand 
their structure. 

Most interestingly, the possible existence of extended critical states in several kinds of 
aperiodic systems, including both quasiperiodic@~0 and non-quasiperiodic one , has been 
discussed in the last few years spurring the interest on the precise nature of CWFs and 
their role in the physics of aperiodic systems. From a rigorous mathematical point of view 
the nature of a state is uniquely determined by the measure of the spectrum to which it 
belongs. In this way, since it has been proven that Fibonacci lattices have purely singular 
continuous energy spectraB, we must conclude that the associated electronic states cannot 
be, strictly speaking, extended in the Bloch's sense. This result holds for other aperiodic 
lattices (Thue-Morse, period doubling) as wellEU, and it may be a general property of the 
spectra of self-similar aperiodie systemS. On the other side, from a physical viewpoint, the 
states can be classified according to their transport properties which, in turn, are determined 
by the spatial distribution of the wave function amplitudes (charge distribution). Thus, 
conducting systems are described by periodic Bloch states, whereas insulating systems are 
described by exponentially decaying wave functions corresponding to localized states. In this 
sense, since the amplitudes of CWFs in a Fibonacci lattice do not tend to zero at infinity 
but are bounded below throughout the systemEl, one may expect their physical behavior to 
be more similar to that corresponding to extended states than to localized ones. 

In this letter we are going to show analytically that a subset of the CWFs belonging to 
general Fibonacci systems are extended from a physical point of view. This result widens the 
notion of extended wave function to include electronic states which are not Bloch functions, 
and it is a relevant first step to clarify the precise manner in which the quasiperiodic order of 
Fibonacci systems influences their transport properties^. We shall begin the investigation 
of this issue by introducing two novel approaches. In the first place, we present a new 
renormalization approach opening, in a natural way, an algebraic formalism which allows 
us to give a detailed analytical account of the transport properties of CWFs for certain 
particular values of the energy. In the second place, we relate the spatial structure of CWFs 
to the topology of the underlying lattice by means of the study of the Fourier spectrum of 
the wave function amplitudes and multifractal techniques. 

The formalism we are going to introduce is based on the transfer matrix technique 
where the solution of the Schrodinger equation is obtained by means of a product of 2 x 2 
matrices. Real-space renormalization group approaches, based on decimation schemes, have 
proved themselves very successful in order to numerically obtain the energy spectrum of 
deterministic aperiodic systemsEJlLa. The convenience for such procedures stems from the 
fact that, by construction, a given transfer matrix relates only three consecutive sites along 
the lattice, so that by decimating the original chain into successively longer blocks we are 
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able to describe the electronic state corresponding to sites more and more farther apart. 
In this context, the key point of our procedure consists of the fact that we renormalize 
the set of transfer matrices instead of the lattice itself. Since these matrices contain all 
the relevant information concerning the dynamics of the electrons, our approach becomes 
specially well suited to describe the characteristic features associated to the long-range 
order of the underlying Fibonacci system for, as we will see below, it preserves the original 
quasiperiodic order of the lattice at any stage of the renormalization process. 

Let us start by considering a general Fibonacci system in which both diagonal and off- 
diagonal terms are present in the HamiltonianEl'Ei 

H = ^2 {V n \n)(n\ + t n , n +i\n)(n + 1| + £ n , n _i|n)(n - 1|}, 

n 

where V n is the on-site energy and t ni „±i are the nearest-neighbor hopping integrals. This 
Hamiltonian can be rewritten in terms of the following matrices 

^w^(^- 1 ), (1) 




where E is the electron energy, a {(3) denote the on-site energies of sites A (B), tAB — tsA 
and tAA are the corresponding hopping integrals and 7 = tAA^AB > 0. Making use of these 
matrices we can translate the atomic sequence ABAAB . . . describing the topological order 
of the Fibonacci lattice to the transfer matrix sequence XYZXWXYZXYZXW . . . describing 
the behavior of electrons moving through it. In spite of its greater apparent complexity, we 
realize that by renormalizing this TMS according to the blocking scheme = XYZ and 
Rb = XW, we get the considerably simplified sequence RaRbR^R^Rb • • • The subscripts 
in the Rs matrices are introduced to emphasize the fact that the renormalized TMS is 
also arranged according to the Fibonacci sequence and, consequently, the topological order 
present in the original lattice is preserved by the renormalization process. Let N = F n be 
the number of lattice sites, where F n is a Fibonacci number obtained from the recursive 
law Fk = Fk-i + Fk-2, with Fx = 1 and Fq = 1. It can then be readily shown that the 
renormalized TMS contains = F n _ 3 matrices R^ and = -F n „ 4 matrices R#. 

We will now use several properties of the Rs matrices to develop our procedure. Firstly, 
they are unimodular (i. e. their determinant is one) for any choice of the system parameters 
and for any value of the electron energy. Secondly, they commute for certain values of the 
energy. In fact, after some algebra we get 



[R A , R B ] - , (2) 




where we have defined the origin of energies in such a way that (3 = —a and tAB = 1, and 

P(E) = ( 7 2 - 1)(E + a)[{E - af - uo 2 } + 2c^ 2 7 2 , (3) 

with uo 2 = (1 + 7)£^g. This commutator considerably simplifies for the two cases mostly 
discussed in the literature, namely the on-site (7 = 1) and transfer (a = 0) models. The 
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expression fl2|) shows that the on-site model is intrinsically non-commutative, for the com- 
mutator vanishes only in the trivial periodic case. On the contrary, there exist three energy 
values for which the R matrices commute in the transfer model, corresponding to E = and 
E = ±u. Most interestingly, since P(E) is a real cubic polynomial in E, there exists at least 
one energy satisfying the relation P(E) = for any realization of the mixed model (i. e. for 
any 7 value). For these energies the condition [Ra,Rb]^= is fulfilled and, making use of 
the Cayley-Hamilton theorem for unimodular matrices^ the global transfer matrix of the 
system, M(N) = R^ A R^ S , can be explicitly evaluated in terms of Chebyshev polynomials 
of the second kind, U m (x), where x is Tr(R)/2 . Alternatively, the required power matrices 
can be evaluated by diagonalizing them to a common basis. From the knowledge of M(N) 
the condition for the considered energy value to be in the spectrum, |Tr[M(iV)]| < 2, can be 
readily checked and, afterwards, relevant parameters describing their transport properties 
can be determined explicitly. In this way, given any arbitrary Fibonacci lattice, we are able 
to obtain a subset of its energy spectrum whose eigenstates can be studied analytically. 

Although our approach is completely general, the solutions of the algebraic equation 
P(E) = are rather involved in most cases. For convenience, model parameters for which 
the algebra considerably simplifies will be discussed henceforth. In Fig. [I] we show the charge 
distribution of two electronic states, for which the R matrices commute, corresponding to 
Fibonacci lattices with N = F 16 = 1597. Figure |l|(a) corresponds to the energy E = 2 and 
lattice parameters 7 = 3 and a — 1. Figure |l](b) shows the result for the energy E = a/1 + 7 
for a transfer model (a = 0) with 7 = 1 + a/2, the so-called silver mean. The overall periodic- 
like behavior of the wave function amplitudes, which we have calculated exactly with the 
aid of our matrix formalism, clearly indicates their extended character. 

In order to discuss this point from a more rigorous physical perspective we will focus our 
attention on the transmission coefficient, t(E). For the transfer model the global transfer 
matrices corresponding to the energies E = ±00 in the interval < 7 < 3 can be expressed, 
after lengthy algebra, in the closed form 



m(n, ±.) = R7= jr,L . (4) 




where p = 2n# — ua = F n s and 2cos# ± = From expression (Q) we get Tr[M(iV)] = 
2 cospO^ 1 and, consequently, we can ensure that the energies E = ±u belong to the spectrum 
in the quasiperiodic limit (N — > 00). Now, we proceed to the calculation of the transmission 
coefficient by embedding the Fibonacci lattice in an infinite periodic arrangement of identical 
atoms connected by hopping integrals t. In this way we obtain 

T(±UJ) = 1+ (5) 

Two important conclusions can be drawn from this expression. In the first place, the trans- 
mission coefficient is always bounded below for any lattice length, which proves the true 
extended character of the related eigenstates. In the second place, we observe that the 
transparency condition r = 1 is obtained, for any value of t, for certain chain sizes satisfying 
pO^ = kir, k = 1,2..., which in turn implies ±00 = =f2 cos(/c7r/p). In this way the trans- 
parent states r = 1 can be classified according to a well defined scheme determined by the 
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integers k and p. Thus the state shown in Fig. [T|(b) corresponds to the choice k = 1 and 
p = 8. 

After having discussed the transport properties of this class of critical wave functions, 
we turn our attention to their spatial structure. To this end we shall consider the example 
shown in Fig. which corresponds to the energy E = V2a, where a = (y/E + l)/2 is the 
so-called golden mean, and model parameters 7 = \fh and a = 0. In Fig. ||](a) we show the 
overall charge distribution through a lattice with N = F 2 \ = 17 711 sites. From this plot we 
notice the existence of two different superimposed structures. In fact, a periodic-like long- 
range oscillation with a typical wavelength of about 18 000 sites is observed to modulate 
a quasiperiodic series of short-range minor fluctuations of the wave function amplitude, 
typically spreading over 122 lattice sites. This qualitative description receives a quantitative 
support from the study of its Fourier transform, as it is shown in Fig. 0(b). In fact, we observe 
two major components in the Fourier spectrum corresponding to the low and high frequency 
regions, respectively. In the low frequency region two relevant features at frequencies V\ ~ 
1.11 x 10 -4 and v<i ~ 8.18 x 10~ 3 are present, in agreement with the short- and long-range 
structures of the charge distribution observed in Fig. 0(a). On the other side, in the high 
frequency region we observe a series of features grouped around frequency positions given by 
successive powers of the inverse golden mean a' - 1 = (V5-l)/2. These feat ures are labelled 
correspondingly in Fig. 0(b) and are characteristic of the quasiperiodic nature of the charge 
distribution of the considered states. By this we mean that we have observed that their 
weights in the Fourier spectrum increase and their positions approach the reported powers 
of a" 1 as the system size increases. 

Finally, to gain further insight into the behavior of the wave function at all length scales 
we have performed a multifractal analysis of the states belonging to the subset E = ±u. 
The amplitude distribution of the electronic states has been characterized by the scaling 
of moments n q {N) of order q, associated to their charge distribution, with the system size 
(for a definition of those moments see, e.g, Ref. |2T]). The multifractal dimension D q is 
determined via the scaling fi q (N) ~ j\f( 1 -i) D <i for q ^ 0. In all cases studied we have found 
that D q — 1, for all q, and for system sizes as large as N = F30 = 1 346 269. Thus the lack 
of multifractality along with the fact that D q equals the spatial dimension proves our claim 
that these states uniformly spread over the whole system. 

In summary, in this letter we prove the existence of a subset of the singular continuous 
spectrum characteristic of Fibonacci systems whose eigenstates are extended in the physical 
sense previously discussed. This we have shown by means of a transfer matrix renormaliza- 
tion technique which allows us to unveil the effects of short-range correlations by grouping 
ABA sites and AB sites into the matrices and R^, respectively. In this sense we can 
properly state that these states are characteristic of the quasiperiodic order of the underlying 
lattice. Interestingly we note that similar results concerning extended states in Thue-Morse 
chains have been recently reported in the literature^. We wish to stress that the algebraic 
approach presented in this work can be extended in a straightforward manner to other kinds 
of aperiodic systems based on substitution sequences, and therefore it can be relevant in or- 
der to attain a unified treatment of physical properties of aperiodic systems. In closing, we 
note that the transparency condition r = 1 obtained for extended Fibonacci states is similar 
to that appearing in random dimer modelsil and this fact suggests that sets of extended 
states may arise in general aperiodic systems. 
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FIGURES 

FIG. 1. Electronic charge distribution in Fibonacci lattices with N = and (a) 7 = 3, a = 1, 

E = to = 2 and (b) 7 = 1 + y/2, a = 0, E = u = \fl + y/2. Insets show finer details of the squared 
wave functions. 

FIG. 2. (a) Electronic charge distribution in a Fibonacci lattice with N = F21, 7 = y/5, 

a = 0, and E = uj = \J 1 + y/b. Inset shows finer details of the squared wave function, (b) The 
corresponding Fourier spectrum. Inset shows the lower frequency region. 
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